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Summary 

Force histories resulting from the impact of compliant 
projectiles were determined experimentally. A long instru- 
mented rod was used as the target, and the impact force was 
calculated directly from the measured strain response. Results 
from a series of tests on several different sized impactors were 
used to define four dimensionless parameters that determine, 
for a specified impactor velocity and size, the amplitude, 
duration, shape, and impulse of the impact force history. 

Introduction 

Impact loading is generally categorized as being in either 
the low-, medium-, or high-velocity regime. The distinction 
between the different velocity regimes is most aptly made not 
in terms of absolute velocities, but rather in terms of the type 
of deformation experienced by the impactor and target during 
the impact event (ref. 1). Low-velocity impact is characterized 
primarily by elastic deformation of impactor and target, and 
by a small amount of highly localized yielding in the material 
immediately surrounding the point of contact. Classical elastic 
contact models (refs. 2 to 6) have been applied successfully 
to a wide variety of low-velocity impact problems, and this 
field is quite well developed. Similarly, high-velocity impact 
has been investigated extensively, particularly during the 
Apollo spacecraft era (refs. 7 to 10), and more recently for 
munitions and space station shielding applications (refs. 11 
to 14). 

The least thoroughly investigated impact regime is the mid- 
range. In this regime, the impact velocity is high enough to 
cause finite deformation of the impactor and/or target, but not 
high enough to generate the extreme amounts of heat and the 
resulting material phase transformations that can occur during 
high-velocity impacts. Comparatively few methods are 
available that can predict with sufficient reliability the transient 


force that results from an impact loading in the mid-velocity 
range. This report contributes toward that goal. 

The objectives of this work are, first, to describe a simple 
experimental technique that can be used to characterize the 
force history resulting from the impact of a highly compliant 
projectile undergoing large nonlinear elastic deformation 
during the impact event, and secondly, to use data from a series 
of these tests as the basis for an empirical force-history model 
that describes the force resulting from such an impact. 

In the first section of this report, a detailed description of 
the experimental apparatus and procedure is given. The results 
of a low velocity “calibration test 1 ’ are compared with finite- 
element calculations to demonstrate the validity of the 
approach. Data from a series of impact tests with silicon rubber 
projectiles are then used to define, for specified impactor 
velocit> and size, four nondimensional parameters that can 
be used to determine the amplitude, duration, shape, and 
impulse of the impact force history. 

Symbols 

D impactor diameter 

E Young's modulus 

F contact force 

F 0 maximum contact force 

/ impulse 

K contact stiffness 

m impactor mass 

n exponent 

R sphere radius 

T contact duration 

t time 

t F time at maximum contact force 


1 



V 

impact velocity 

Vo 

reference velocity 

w 

diagonal matrix of weighting factors 

w,- 

curve-fit weighting factor 

X 

vector of abscissa measurements 

y 

vector of ordinate measurements 

a 

relative indentation 

V 

Poisson’s ratio 

n,-n 4 

nondimensional parameters 

7 T 

3.14 

P 

mass density 


Superscripts: 
a,0,y,6 exponents 

Apparatus and Procedure 

Longitudinal stress waves propagate nondispersiveiy in a 
uniform thin bar. An instrumented bar is therefore well suited 
to determining the force history due to longitudinal impact 
(refs. 5 and 15). 

The experimental apparatus is shown in figure 1 . The target 
was a uniform aluminum 6061-T6 bar 91 cm in length, with 
a 1.27-cm-diam circular cross section. Two electrical 
resistance foil strain gages (type EA-1 3-062 AQ-350) were 
mounted in series at diametrically opposing points at the 


midpoint of the bar. This gage arrangement was chosen so 
as to eliminate the effect of any small flexural waves that would 
be generated by a slightly eccentric longitudinal impact. The 
bar was suspended with lightweight 24-gage wire at two 
support points located at 30 and 60 cm from either end. This 
arrangement allowed the bar to swing freely in the longitudinal 
direction after the impact occurred, and resulted in a true free- 
free boundary condition for the bar. A 2.54-cm-diam 
aluminum end cap was installed on the proximal end of the 
bar, as shown in figure 2. The purpose of the end cap was 
to transfer the compressive force from the compliant impactor 
to the bar as the impactor deformed and flattened out to a 
greater diameter than that of the bar. 

The impactors were fired at the bar from an air gun. The 
impact velocity was determined from the measured transit time 
between two photoelectric diodes which were placed 25 and 
50 cm from the exit end of the gun. This sensor arrangement 
may slightly underestimate the actual velocity of the projectile 
when it impacts the target. Silicon rubber balls of two different 
diameters were used as the impactors. Mechanical properties 
of the impactors are given in reference 16 and summarized 
in table I. As the compressive strain pulse generated by the 
longitudinal impact passed the strain gages, the change in 
voltage output was amplified by the preamplifier (shown in 
fig. 1), temporarily stored in the waveform recorder, and 
displayed on the oscilloscope. Permanent copies of the 
recorded signal were made on the plotter, and manually 
digitized records of the data were then stored in a computer 
for later analysis and plotting. 


ELECTRONIC TIMER 



Figure 1 .—Longitudinal bar experimental setup. 
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Figure 2.— Bar instrumentation and fixturing. 


F = Ka n 


(i) 


where 

F = contact force between spheres 
a = relative indentation between spheres (u { - u 2 ) 
n = 1.5 

and 



where 

R ( = radii of spheres 


and E r v { are the respective elastic constants. A special case 
of interest here occurs when the target is flat (R 2 = °°) in 
which case equation (2) simplifies to 




( 3 ) 


TABLE I. -MECHANICAL PROPERTIES 
OF SILICON RUBBER IMPACTORS 

[From ref. 16.] 


Density, gm/cc 1.03 

Hardness, Shore A 65 ±5 

Tensile strength, MPa 5.2 

Maximum elongation, % 110 


Calibration Test: Hertzian Impact 

Before proceeding to the complicated problem of the high- 
speed impact of the compliant projectiles, a simpler case was 
considered. In order to validate the experimental method, a 
preliminary test and corresponding analysis of the impact of 
a steel ball on the aluminum bar were performed. During a 
low-speed impact, the deformation should be mostly elastic, 
and thus the contact behavior between the steel and aluminum 
should be well described by the Hertz model, which is 
summarized briefly here. 

Hertz derived the force-indentation relation to describe the 
elastic contact behavior between two spherical bodies. Hertz’s 
contact law (refs. 2 and 3) is commonly used in both static 
and dynamic applications to determine the force arising from 
elastic contact. To briefly summarize Hertz’s contact law, we 


Figure 3 shows the strain history measured at a single gage 
location on the bar during a low-speed, elastic impact, at 
1 .6 m/s. The incident compressive pulse (100 < t < 150 fis) 
generated by the impact reflects from the free end of the bar 
as a tensile pulse (275 < t < 325 ps) and propagates back to 
the strain gage. For comparison, a finite-element impact 
analysis of this calibration test was performed. The bar was 
modeled with a series of four-degree-of-freedom rod elements 
(ref. 17) which used the longitudinal displacements (u) and 
the corresponding strains ( du/dx ) at both nodal points as the 
degrees of freedom. The Hertz contact law (eq. (1)) was 
incorporated into the finite-element program to define the 
force-displacement relationship between the impactor and the 
bar. Newmark’s implicit method of direct-time integration 
(ref. 18) was used to solve the equations of motion for the 
impactor and the bar simultaneously. The impactor was 
represented by a lumped mass with an initial velocity toward 
the bar. 

The experimental results in figure 4 compare reasonably well 
with the finite-element analysis using a 39-element (80-degrees- 
of-freedom) model of the bar. The contact force history is 
inferred from this strain measurement by multiplying the strain 
by the axial stiffness EA of the aluminum bar and by assuming 
that the initial pulse travels undistorted (nondispersively) 
througli the bar. The impact force determined in this manner 
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Figure 3.— Strain history at midpoint of bar from 1.6 m/s impact of 
1.59-cm-diam steel ball. 



TIME, t, ms 

Figure 4. — Force history from 1.6 m/s impact of 1.59-cm-diam steel ball 
on aluminum bar. 


from the strain in figure 3 is compared with that calculated 
during the finite-element analysis in figure 4. The bar translates 
longitudinally during contact with the projectile, so the force 
history is not quite symmetrical about the maximum force, 
as would be expected if a more rigid target were used. 
Apparently, the finite-element model accurately calculates the 
strain history due to a Hertzian impact. In addition, the 
measured strain response provides a reasonable estimate of 
the impact force. 


Longitudinal Impact With Soft Projectile 

The same bar that was used in the calibration test was then 
used to measure the strain history resulting from the 
longitudinal impact of a 1 .27-cm-diam silicon rubber ball. The 
force history will then be calculated from the strain, and from 
here on will be referred to as the “measured force.” 

Figure 5 shows a typical measured impact force versus time 
behavior for the impact of a 1 .27-cm-diam silicon rubber ball. 
Unlike the more nearly symmetrical Hertzian case, the shape 
of the curve is noticeably skewed. The force reaches its peak 
very early in the contact interval and then tapers off slowly 
before the impactor loses contact with the target. The Hertzian 
curve in figure 4 could be accurately approximated by a simple 
half-sine wave of amplitude F 0 and duration T. The non- 
Hertzian behavior of the soft impactor, on the other hand, is 
non-symmetrical; in addition to F () and T, it requires an 
additional parameter, t Fo (the time at which the peak force 
occurs), to characterize its variation with time. Figure 5 
identifies the three parameters on the measured curve. 

A series of tests was performed to determine how the force 
history in figure 5 varied with impact velocity. Figures 6 to 
8 show the variation in shape and amplitude of the force history 
resulting from impact of the 1 .27-cm-diam balls at velocities 
ranging from 25 to 150 m/s. Figure 6 shows that the amplitude 
of the force varies in proportion to F 2 , which is different 
from the nearly linear variation with velocity predicted by the 
classical elastic impact models (ref. 4). Figure 7 shows that 
the contact time varies inversely with the impact velocity, also 
in contrast to the elastic models, which predict that the two 
are independent. In plotting these data, the “effective contact 
time” T is taken as the duration of the large initial portion 
of the force history pulse, during which the majority of the 
impulse is transferred to the target. This approximation is 
shown in figure 5, and has the effect of ignoring the long 
trailing “tail” of the force-history curve. 

In Figure 8, the time t F at which the maximum impact 
force occurs is plotted over the velocity range tested. As the 
impact velocity is increased, the maximum force is reached 



Figure 5.— Force history measured from 88.9 m/s impact of I.27-cm-diam 
silicon rubber impactor. 
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earlier in the contact interval until it remains approximately 
constant at velocities above 65 m/s. 

The impulse transferred to the target during contact with 
the impactor is 

/= |V(r) dt (4) 

where F(t) is the time history of the impact force and T is the 
contact duration. Assuming that the force history in figure 5 
can be approximated by two sinusoidal waves as shown, we 
have 


7T 

F n sin — 
‘ 2 



for 0 < / < t F 


F(t) = 


7T 

cos — 
2 



for t F < t < T 


(5) 



Figure 7.— Variation of contact time with impact velocity. 


0 for t > T 

which, after substitution into equation (4), give 

/= -FJ (6) 

7T 

The impulse calculated by using equation (6) is shown in 
figure 9 to vary linearly with impact velocity. 

A second series of tests was conducted with 0.95-cm-diam 
impactors of the same silicon rubber material. The analogous 
results are overlayed on figures 6 to 9 for comparison. The 
trends observed in the earlier series of tests apparently do not 
depend on the impactor size, but instead are characteristic of 
the material. All measured data show the same trends. The 
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Figure 9. — Measured impulse versus impact velocity. 
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numerical differences in the measured impact force histories 
are due to the smaller mass of the second impactor. 

Dimensional Analysis 

An empirical model is now developed to describe the 
observed variation of impact force history with impactor size 
and velocity. The development and notation follow that of 
Buckingham’s “Pi Theorem” (ref. 19). 

Specifically, we wish to describe the data in figure 6 , which 
show the variation of maximum contact force F 0 with impact 
velocity V and impactor size, D. We therefore define a 
parameter II j, as: 

n, = FS p 8 v D 6 (7) 

where p is the mass density of the impactor, and the exponents 
are chosen such that II j is dimensionless. The exponents are 
determined by rewriting equation (7) in dimensional form, that 
is. 




shown in figure 10 is calculated by using a least-squares curve- 
fit (ref. 20 ): 


By equating exponents of M on the left side of equation ( 8 ) 
with those on the right, we have 

0 = a + (3 (9) 

and by doing the same for L, we have 

O = a-30 + 7-h<5 (10) 

and similarly for T: 

0 = 2a — 7 ( 11 ) 

In solving equations (9) to (11) simultaneously for the four 
exponents, we may choose a = 1 for convenience, giving 


m — 


x T Wy 

x T Wx 


(15) 


where, if n separate measurements have been taken, x and y 
are (n X 1) vectors of the ordinate and abscissa values, 
respectively, and 


W 


W] 

w 2 

w 3 


(16) 


[a, jS, 7 , 6 ] = [1, -1, -2, -21 (12) 

or, using the definition of II j in equation (7), we have 


n, 


Fq 

pV 2 D 2 


(13) 


as the dimensionless parameter. The data in figure 6 have 
therefore been replotted in figure 10 by using the denominator 
( P V 2 D 2 ) in equation (13) as the X axis. A linear curve-fit is 
used to calculate the slope of the line that best fits the measured 
data. The slope m of the line 

y = mx ( 14) 


is an ( n x n) diagonal matrix of weighting factors. Figure 10 
shows that the variability in the measured force increases as 
the value of pV 2 D 2 is increased. The weighting function for 
the linear curve-fit shown in figure 10 is therefore taken as 

2 

i=\,n (17) 



where V x is the impact velocity for test “i” and V () is a 
constant reference velocity arbitrarily chosen to make w, 
dimensionless. As a result, the data for the higher impact 
velocities are weighted less heavily that those at the lower 
velocities, giving the curve-fit a closer correlation in the low- 
velocity range, where the linear approximation is more valid. 
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According to equation (13), the slope of the line in figure 10 
can be interpreted as an estimate of Ilj. 

The data plotted in figure 7 show the variation of contact 
time T with impactor size and velocity. Following the 
procedure outlined above, we define for this case a second 
dimensionless parameter n 2 , such that 


where t F has units of time, and the values of the exponents 
are to be determined. Repeating the above process for n 3 
gives 


n 3 



( 22 ) 


n 2 = T a pp V y D b (18) 

where Thas units of time, and the remaining variables have 
the same definitions as previously. By repeating the process 
used above to determine the dimensionless form, we have 


n 2 


TV 

D 


(19) 


The data in figure 7 are replotted in figure 1 1 by using the 
ratio D/V as the X axis. The data vary linearly at the higher 
impact velocities (lower values of D/ V) so the least-squares 
weighting function is taken as 


so the data in figure 8 are replotted as a function of D/V in 
figure 12. The data vary linearly throughout the range of 
impact velocities tested, so no weighting is used in the linear 
curve-fit The slope of the line in figure 12 is an estimate of 
n 3 for the given impact conditions. 

Finally, a fourth dimensionless constant is used to describe 
the variation of impulse / with impactor size and velocity. We 
define in this case 


n 4 = /“ p' 3 v y d 6 


(23) 


w,= 



where impulse / has units of (ML/T) and the exponents are 
(20) to be determined. The procedure used previously gives 


in figure 1 1 . This allows the linear curve- fit to most accurately 
approximate the data in the region where the linear assumption 
is most valid. From equation (20), the slope of the line in 
figure 11 is an estimate of II 2 . 

The data plotted in figure 8 show how the time of maximum 
force, t F , varies with impactor size and velocity. To describe 
this variable, we define a third dimensionless parameter, n 3 , 
such that 

U 3 = t a F yV y D b ( 21 ) 


n 4 = 


i 

pVD 3 


(24) 


The impulse data in figure 9 are therefore replotted in figure 1 3 
with the value of pKD 3 as the X axis. The slope of the line 
is an estimate of the parameter n 4 . Because the impulse is 
calculated directly from the measured values of F 0 and T 
according to equation (6), it follows that n 4 is not inde- 
pendent of the previous three dimensionless parameters. 




Figure 12.— Least-squares curve-fit for n v 
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Figure 13.— Least-squares curve-fit for n 4 . 

Indeed, the substitution of equations (6) to (8) into equation (10) 
shows that 


7 T 

as can be verified with the values given in table II. 

The measurements of F () , T, t Fo and /, (shown graphically 
in figures 6 to 9) are, of course, specific to the particular 
impactor /target combination used in these tests. Therefore, 
the values of the four Pi parameters, summarized in table II, 
will be different from those given here when equations (7) to 
(10) are applied to different impactor/target combinations. 
Nonetheless, for impact scenarios involving a large nonlinear 
elastic deformation of the projectile, the nondimensional 
modeling approach outlined in the previous section can be used 
to extrapolate or “scale” the impact force measured from a 
single test over a range of impactor velocities and sizes. 

Conclusions 

A simple experiment to measure the force history resulting 
from the ballistic impact of a compliant projectile is described. 


TABLE II.— DIMENSIONAL ANALYSIS FOR 
IMPACT FORCE HISTORIES 


[Least- squares curve-fit data for figures 10 to 13.] 


Parameter 

Dimensionless 

ratio 

Empirical 

value 

Weighting 

factor, 

w i 

Reference 

figure 

number 

n, 

F„ 

1.096 

/v\ 2 

10 


pV^D? 


\vj 



TV 


fy, V 


n 2 

— 

0.946 

( — ) 

11 


D 


w 



h/ 




n. 

— 

.231 

1.0 

12 


D 





l 




n 4 

— 

.672 

1.0 

13 


pVD 3 

n 4 . 




The results show that for the compliant rubber impactor used 
here, the velocity dependence of the force history is sig- 
nificantly different from that which would be predicted by 
classical impact mechanics. This is due primarily to the large 
nonlinear elastic deformation of the impactor. From the test 
results, four nondimensional parameters are determined that 
can now be used to estimate the force history and impulse 
transferred by a soft impactor to the target for a given set of 
impact conditions. 


Lewis Research Center 

National Aeronautics and Space Administration 
Cleveland, Ohio, June 13, 1990 
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